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ABSTRACT 

We use new Herschel multi-band imaging of the Andromeda galaxy to analyze how dust 
heating occurs in the central regions of galaxy spheroids that are essentially devoid of young 
stars. We construct a dust temperature map of M31 through fitting modified blackbody SEDs 
to the Herschel data, and find that the temperature within 2 kpc rises strongly from the mean 
value in the disk of 17 ± 1 K to ~ 35 K at the centre. UV to near-IR imaging of the central 
few kpc shows directly the absence of young stellar populations, delineates the radial profile 
of the stellar density, and demonstrates that even the near-UV dust extinction is optically thin 
in M31's bulge. This allows the direct calculation of the stellar radiation heating in the bulge, 
U*(r), as a function of radius. The increasing temperature profile in the centre matches that 
expected from the stellar heating, i.e. that the dust heating and cooling rates track each other 
over nearly two orders of magnitude in U*. The modelled dust heating is in excess of the 
observed dust temperatures, suggesting that it is more than sufficient to explain the observed 
IR emission. Together with the wavelength dependent absorption cross section of the dust, 
this demonstrates directly that it is the optical, not UV, radiation that sets the heating rate. 
This analysis shows that neither young stellar populations nor stellar near-UV radiation are 
necessary to heat dust to warm temperatures in galaxy spheroids. Rather, it is the high densities 
of Gyr-old stellar populations that provide a sufficiently strong diffuse radiation field to heat 
the dust. To the extent which these results pertain to the tenuous dust found in the centres of 
early-type galaxies remains yet to be explored. 
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1 INTRODUCTION 

As the nearest, massive galaxy, Andromeda (M31, NGC224) has 
offered a unique insight into the properties of galaxies. It provides 
the perfect stepping stone between the well-resolved interstellar 
medium (ISM) and stellar populations within our own Galaxy, and 
the integrated properties of more distant galaxies. Due to its prox- 
imity (~ 780kpc), Andromeda offers more than a resolved exam- 
ple of an early-type spiral galaxy, but can also be used to explore 
the interaction of the ISM and stars in early-type galaxies (ETGs) 
through its large, visually dominant bulge. 

Andromeda's bulge dominates the stellar luminosity and 
mass within the cen tral 1.5 kpc (K c ff(bulge) ~ 0.5 - 1.1 kpc; 
ICourteau et al.ll2oTlh . with the disk dominating beyond this. The 
bulge contributes ~30% of the total stellar mass and luminosity in 
M3 1 , and it is clearly the highest surface brightness feature at UV - 
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NIR wavelengths dGeehan et alj200fj ; [Courteau et alj2oTTh . While 
the integrated colours and l uminosity of M3 1 may place the galaxy 
in the "green valley" (M utch et alj201 lb . t he optical colours of the 
bulge are red, with a B - V ~ 0.9 to 1.0 dWalterbos & Kennicuttl 
Il987t) . placing it securely in the rang e of the red-sequence, where 
most early-type galaxies are found. lOke & Sandagel (1968) even 
used the spectral energy distribution of the centre of M3 1 as a rep- 
resentative for the average giant elliptical galaxy when determining 
^-corrections. 

This "early-type" nature of the bulge is supported by the old 
mean stellar age determined for the central region of M31. Re- 
solved star colour-magnitude diagrams created by ground-based, 
adaptive-optics, NIR imag ing reveal a popul a tion dominated b y 
stars greater than 6 Gyr (IDavidge et alj|2005l: blsen et al.ll2006h . 
while single-stellar population fi ts to absorpt i on-lin e indices from 
slit-spectroscopy of the bulge ( Sagli a et all |2010|) find that the 
bulge of M31 is uniformly old (> 12 Gyr, excluding the cen- 
tral arcsecs). High-resolution, individual star photometry from 
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the Pan-Chromatic H ubble Andromeda Treasury survey (PHAT; 
iDalcanton et alj|2012l) find no population of young-stellar sources, 
with t he UV-light dominated b y evolved stars such as post-AGB 
stars feosenfield et alJlin prep.l) . 

Given this lack of young stars and star-formation, it is unsur- 
prising that the bulge of M3 1 is also ex tremely gas poor. L ittle to no 
CO is detected in the CO(l-O) map of lNieten etal] d2006h down to 
low surface brightnesses. CO is detected in the central part of M31 
when deeper observations of small high attenuation r egions are 
made (Mel chioret all 2000; Me lchior & Combesll2oT ]]), however 
the attenuation in these regions is still relatively low (Ab < 0.3) and 
the co vering fraction of these regions is very small dMelchior et al.l 
2000), meaning diffuse, low attenuation dust (and presumably dif- 
fuse gas) is more c haracteristic of the bulge. Compounding this is 
a low CO-H 2 ratio dLerov et"aill201 lb . and a low total HI column 
jBraun et alj2009h giving a total cool gas mass in the centre of only 
a few 10 6 M o (~ 0.02 % of the bulge s tellar mass), likely dominat ed 
by molecular gas (e.g. lli et al .1120091 ; iMelchior & Combesll201l|) . 

This low gas fraction makes the M31 bulge more gas-poor 
than many ETGs, as recent work with the ATLAS 3 D sam ple 
dCappellari et alj|201 lh has demonstrated. lYoung et al"] (1201 lh ob- 
served CO(1-0) in ~22% of the sample of nearby early- type galax- 
ies, giving corresponding gas masses of MQ&i) > 10 7 M o (for a 
sample with a median stellar mass of M t = 3 x 10 10 M o ). Similarly, 
studies of nearby field ellipticals and lenticulars have found atomic 
gas in a large fraction of these (~ 70%), wit h this H 1 emission gen - 
erally associated with ionized gas emission (Morganti et al. 2006). 
Interestingly, the molecular gas in early-type galaxies is not always 
associated with star formation. Optical colours and emission-line 
ratios indicate other forms of heating in the ISM o f some early- 
types , and only low levels, if any, of star formation dCrocker et al.l 

Eou|). 

Given the very low level of star formation, low attenuation, 
and small amount of gas at the centre of M31, little dust emission 
was expected in the central kiloparsecs of M 3 1 . Yet when this re - 
gion was examined in the far-IR with IRAS dHabingetalJI 19841) . 
ISO dHaas et alj [r9981, and Spitzer dGordon et al.ll2006l) . emission 
at wavelengths greater than 60 pm was clearly seen, indicating the 
presence of warm dust. In n earby ETGs, low-l evels of warm dust 
are also being detected, with lSmith et al.l J201 lh finding dust emis- 
sion in 24% of ellipticals and 62% of SOs in the Herschel Reference 
survey. While M ias JM, is far lower in spheroids than in disks, the 
dust appears to be warmer on average than that found in later- type 
galaxies. A similar result was fo und using Hersche l in n earby early- 
type spirals similar to M3 1 by Engelbracht et al.l d2010l) . with the 
mean dust temperature of the bulges consistently hotter than the 
disks in these galaxies. Relatively h igher dust temperatures were 
also found by Rowla nds et alj d201lh in the Herschel- ATLAS sur- 
vey of more distant ellipticals. Hence dust, when detected, tends to 
be in a warmer state in spheroids than that found across the disks 
of later-type galaxies. 

Given the high density of stars in bulges, stellar heating is the 
likeliest explanation for the observed warmer dust temperatures in 
ETGs. Yet the old stellar ages found in ETGs, and especially in the 
centre of M3 1, argue against the standard view that the I R luminos- 
ity an d warm dust is a direct tracer of star formation (e.g. iKennicuttl 
1998). Based on the IRAS observations and the extremely weak 
UV observed using t h e Astronomica l Neth erlands Satellite (ANS; 
IColemanetal]|l980h . lHabing et al] dl984h put forward the argu- 
ment that it is the high density of late-type giant stars that provide 
the strong enough radiation field to heat the dust. Based on Her- 
schel PACS and SPIRE maps of M3 1 with unprecedented resolu- 
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M31 nucleus position b 00* 42™ 44.35" 

(J2000) +41°16'08.60' 

Position angle of major axis 37.7° 

Inclination 75° 

Distance 780 ± 40 kpc d 

Distance modulus 24.46 



a Based on N ED data and references where given. 
b tevans et al.ll2010l) . and verified in Spitzer IRAC 

3.6/im image. 

° IStanek & GarnavictJ fl993) ; lRich et a?]d2005l) . see 

also NED for other determinations. 
d V = 111 ± 12pc along major axis 

tion (Krause et al., in prep.), we follow lHabing et alj d 19841) and use 
this high resolution to demonstrate how and by what stellar popu- 
lations dust is heated in the bulge of this galaxy, and, by proxy, in 
the spheroids of other inactive early-type galaxies. 

We briefly review the Herschel data and reduction in section|2] 
and discuss the FIR geometry and integrated properties of the M31 
bulge in section|3] In section|4]we determine and discuss the mean 
dust temperature across M3 1 and in the centre, and determine the 
heating mechanism for the rise in dust temperature in the centre. 
We finish with the discussion and summary in sections [5] and [6] 
The global properties of M3 1 we assume throughout this paper are 
listed in TableQ] 



2 DATA AND REDUCTION 

The M31 Herschel data and their reduction are described and dis- 
cussed in detail in Krause et al. (in prep.), so we only briefly review 
the data here. 

M3 1 was imaged in all 6 Herschel photometric bands (PACS 
70, 100, and 160 pm, and SPIRE 250, 350, and 500 pm) in slow 
parallel mode for a total time of ~ 24 hours. The images extend for 
~ 3° x 1°, centred slightly off-nucleus (at the PA. given in TableQ}, 
covering the main stellar disk of Andromeda, including the lOkpc 
ring. 

All images were reduced to level one using HIPE version 6.0, 
and then used SCANAMORPHOS vl2.0 dRousse<l20~il) to pro- 
duce the final images. As HIPE v6.0 was used, we converted the 
PACS images from flight model (FM), 5 to FM, 6 by dividing 
by the factors listed in PACS Photometer Point Source Flux Cal- 
ibration Report vl.0, which are 1.119 (70pm), 1.151 (100pm), 
and 1.174 (160 pm) for the three PACS photometer bands, respec- 
tively. All images were converted to MJysr -1 . Based on SPIRE 
Observers' manual v2.4E we used beam areas of 423n" (250 pm), 
75 In" (350 fim), and 1587n" (500 pm) for the three SPIRE bands, 
respectively, to convert from Jy beam -1 to MJy sr~' . 

The mean FWHM of the PACS and SPIRE Point Response 
Functions/beams are; ~ 5.6" (70pm), ~ 6.8" (100pm), ~ 11.4" 
(160pm), 18.2" (250pm), 24.9" (350pm), and 36.3" (500pm), re- 
spectively, for the 20"/s scans used here (for full details of the cor- 
rect point spread functions and beams, see the PACS0 and SPIRE 1 
Observer manuals). The actual images use pixels sizes of 1" for 
all PACS bands and 6", 10", and 14", for the respective SPIRE 

1 http : //herschel . esac . esa . int/Docs/SPIRE/html/spire_om . html 

2 http : //herschel . esac . esa . int/Docs/PACS/html/pacs_om . html 
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bands. When comparing bands, we convert all bands to the lowest 
resolution and pixel size using the convolution kernels provided by 
lAniano etailfcOllh . 

For all bands we have measured the noise and background as 
discussed in Krause et al. (in prep.). The 'background' in the im- 
ages is a combination of both Galactic foreground and background 
galaxies, and is significant only in the SPIRE images. For simplic- 
ity, a uniform background in each band was assumed and subtracted 
from the image. 

The measured noise for each band in the original images 
is; 3.17 MJysr 1 (70jian), 3.23 MJysr 1 (100 pm), 2.29 MJysr 1 
(160pm), 0.68 MJysr 1 (250pm), 0.46 MJysr 1 (350pm), and 
0.20 MJy sr 1 (500 pm). When convolved and binned to lower res- 
olutions, the noise per resolution element in each band is obviously 
reduced, with the effective resolution driven by the longest wave- 
length (e.g. 36.3" FWHM and 14" pixels at the 500 pm resolution). 



3 THE BULGE OF M31: STARS, DUST, AND GAS 
3.1 The Far-IR emission in Andromeda 

In Figure Q] we show a RGB image of the full disk of Andromeda 
using the PACS 70 pm, PACS 100 pm, and SPIRE 250 pm bands. 
One of the things that stands out in this image, apart from the dusty 
ring, is the blue centre of the IR map, indicating a distinct rise in 
dust temperature towards the centre of M3 1 . 

The central kiloparsec bulge region (hereafter referring to a 
circular aperture of radius 4.405 ' ( 1 kpc), centred on the stellar peak 
defined in Table QJ is relatively weak in the SPIRE 250 pm band 
compared to the extended disk and star forming ring, but shows a 
relative excess in the PACS 70 pm band. The far-IR emission in this 
central region appears more circular and shows a spiral like pattern 
within. 



3.2 The correlation of dust and gas 

The spiral pattern can be seen in finer detail with Figure [2] where 
we show the PACS 70pm and SPIRE 250pm c ontours overlaid 
on the Ha image of the c entral kiloparse c s from iDevereux et al.l 
dl994l) . As demonstrated bv lLi et aTI (l2009) using Spitzer data (par- 
ticularly their Figures 6 and 8), the IR emission follows well the 
morphology of the Ha emission. The correspondence is not linear 
between the Ho- and the two Herschel bands, with offsets between 
the Ha peaks and IR peaks, but in general the same barred spiral 
pattern is seen. The two Herschel bands match very well with each 
other within the inner 0.5 kpc (2.2') radius, but outside this the cor- 
respondence becomes weaker as the dust becomes cooler, and the 
PACS 70 pm emission merges with the noise. To the NW and SE 
of the image in Figure |2p) we can see the contribution of the disk 
appearing at the SPIRE 250 pm wavelengths. 

The correspondence of the dust and ionized gas emission to- 
gether suggest that the dust is associated with the gas a nd also likely 
to be in the same lower inc lination thin disk spiral ( Jaco bv et al.l 
1 19851 : ICiardullo et"aT]| 19881) . The offset between the peaks in the 
Ha and the IR emission indicates the presence of cooler, denser gas 
associated with the dust, that is still too low in column density to be 
observed in H i or CO, with the weak and non-dete ctions in CO of 
severa l hig her attenuation regions in th e centre by iMelchior et alj 
J200d) and Melchior & Combe] Hjl 1) supporting this. As the gas 
emission indicates that most of the gas is in a diffuse state, and at- 
tenuation maps show an, on average, extremely low attenuation, it 
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Figure 3. The spectral energy distribution (SED) of the central kpc region 
of Andromeda, measured within a circle of lkpc (4.405') radius from the 
centre. The fluxes (GALEX FUV and NUV, SDSS u, g, r, i, and z, Spitzer 
IRAC and MIPS24, and all six Herschel bands) within this aperture are 
marked by the red points, with the error bars including calibration, sky and 
noise uncertainties. The black curve shows the best fi t model SED from the 
SED fitting code MAGPHYS Ida Cunha et"al]|2008h . while the blue curve 
shows the corresponding unattenuated stellar SED from the same model. 
Below the model SED, the residuals between the observed data and models 
are shown, with |log(L„|, s /L mo( j)| < 0. 1 dex for all wavelengths except for 
SPIRE500pm, where the observed flux is under predicted by ~30%. 



is clear that the dust is likely to be optically thin to a significant 
amount of the radiation from the bulge stars. 

3.3 The spectral energy distribution of the bulge 

The low optical depth in the centre of M3 1 is also made clear by 
the full UV-IR SED in the central kilopars ec, which is illustrate d 
in Figure[3] drawing on data from GALEX (Gil d e Paz et al.i 2007). 
SDSS ( Aih araetalhoilh , Spitzer IRAC tearmbv et alj|2006f) and 
MIPS l IGordon et aiX2 006). and now Herschel. Overlaid on the ob- 
served data is the bes t fit from the SED fitting code MAGPHYS^ 
dda Cunha et al. I l2008l) (black curve) and its implied unattenuated 
stellar emission (blue curve). For this model fit we c r eated a new 
stellar library using the original iBruzual & Charlotl (120031) code 

da Cunha et alT(l2008l) 



and the same star formation histories as in 
but excluding all models which were f ormed less than 6 Gyrs ago 
(r„ < 6Gyr, see Ida Cunha et all d2008l) , specifically §3.1.1), based 
on previous estimates for the mean stellar age of the bulge o f M3 1 
(see section[T] and e.g. iDavidge et al.l2005l : [Olsen et alj200rj) . The 
MAGPHYS model fit returns a stellar mass within this lkpc ra- 
dius of log(M*/M ) = 10.01 + 0.01, and an unattenuated stellar 
luminosity of ~ 10 99 L G . The dust luminosity of the M31 centre is 
~ 10 8 3 L Q , with MAGPHYS determining a total dust mass of only 
log(M dusl /M ) = 5.17 ±0.05. 

As stated in the introduction, the 10 10 M o stellar mass of the 
bulge makes it a significant frac tion (~ 30%) of the total stellar 
mass of M31 dGeehan et al.l2006l) . The optical-UV colours clearly 
reveal the old stellar age of the bulge, with NUV-r « 5.0. These 
colours lead to an extremely low estimated specific star formation 
rate from the model, with sSFR< O.OlGyr 1 (SFR < 10~ 2 M o ), 
in agreement with the lack of young stars observed in the bulge 



' www.iap.fr/magphys 
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Figure 1. Far-infrared image of Andromeda using the Herschel bands showing PACS 70 pm (blue), PACS 100 /im (green), and SPIRE 250 /im (red), all 
convolved to SPIRE 250/im resolution (from Krause et al., in prep.). All three bands have the same square root scaling from 10 MJy sr~' to 150 MJy sr~'. The 
angular scale of the image is shown by the 2.3 kpc (10') bar in the lower left. 
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Figure 2. Continuum-subtracted Ha+ [Nn] narrow band image of the central kiloparsecs of M3 1 (grey scale) from Devereux et al 1 fl99i overlaid with contours 
from PACS 70 pm (left) and SPIRE 250pm (right) imaging. In both images the Ha has a square root stretch, and the centre is marked by a cross of 1 kpc 
length along both the major (4.4') and minor axis (1.14'), assuming the inclination of M31 given in Table[T] For the PACS 70 pm contours, the levels are 15, 
30, 70, and 100 MJy sr~', while for the SPIRE 250 pm contours, the levels are 10, 15, 20, and 25 MJy sr~', each indicated by an increasing thickness. 



of M31 . dDavidge et aill2005l : lolsen et alJl200Q : iRosenfield et al.l 
lin prep.1) . 

It should be noted that the values here for the dust mass and 
luminosity in the central kiloparsecs are likely upper limits for the 
dust in the bulge, as even at these small radii, some fraction of the 
IR emission still arises from the disk in the foreground and back- 
ground due to the high inclination of Andromeda. This contribu- 
tion of the disk infrared emission may also explain the weak UV 
(< 0.25 pm) excess observed in the unattenuated spectrum of the 
best fit MAGPHYS model relative to the observed SED. To explain 
the observed IR SED, the MAGPHYS model finds a small com- 
ponent of buried young stars necessary. This component is small 
(as demonstrated by the low sSFR), and cannot be associated with 



the bulge (as the high-resolution observations have shown), but still 
contributes significantly to the unattenuated UV. 

The total attenuation in the diffuse ISM is extremely low, as 
can be seen in the differences between the black and blue curves, 
with the model finding an attenuation in the diffuse ISM of only 
Tv ~ 0.03. The low value of attenuation measured from the SED 
match es that observed in B-band attenuation maps (Melchi or et al.l 
2000), and is a result of a very low dust column, with the SED fit re- 
turning an average dust column density of only ~ 4.5x 10 4 M o kpc~ 2 
in the central region. 

In contrast to the stellar mass, the bulge contributes little to 
the total dust mass in M31. At larger radii, emission from the disk 
of M31 along the minor axis overlaps with (and dominates over) 
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the dust emission from the bulge. Fits to the integrated IR spectral 
energy distribution (SED) of the Andromeda galaxy ( within a 21.5 
kpc ap ertu re, Krause et al . , in pr ep.), using both the iDraine & Lil 
J2007I) and lda Cunha et all {2008) models, reveal that the dust mass 
in the inner 2kpc only contributes ~0.5% to the total dust mass of 
M3 1 . However, the central 1 kpc contributes a ten-times larger frac- 
tion (i.e. ~ 5%) to the total IR luminosity of M31, because of the 
relatively warmer dust emission. We can quantify the temperature 
of this dust in the bulge and in the rest of M3 1 by fitting simple 
modified blackbodies. 



4 THE DUST HEATING IN THE M31 BULGE 
4.1 Simple modified blackbodies 

We now convert the spatially resolved thermal-IR SEDs into dust 
temperatures, describing the emission locally with a single dust 
temperature, T d . While a simplification, this method returns a rea- 
sonable estimate for the mean, luminosity-weighted temperature of 
the dust. 

We do this by using a simple modified blackbody; 

S v = N 6 K v0 B v (J i )(^j , (1) 

where the dust surface brightness, S v , is proportional to the Planck 
function for the given dust temperature, i? v ,(T d ), modified by the 
dust emissivity, which is assumed to be a power-law function of 
frequency, Kyo (v/vof- At a given dust temperature this is then sim- 
ply multiplied by the dust column density, 2V d , to yield the surface 
brightness. 

The power-law fun ction is a reasonab l e approximat ion for the 
dust emissivity (see e.g. Hildebrand 19831: lDraindl2003h . with the 
alactic diffuse ISM well fit by an emissivity index of /? = 2 
Draine & Ledl 1984b . While a mean temperature for the dust can 
be determined, in reality there will be a range of temperatures along 
the lines-of-sight due to both a distribution of gra in sizes and heat- 
ing radiation fields (see, e.g. IDraine & Lill2007h . This simplified 
modelling foremost reflects a physically-motivated conversion of 
the FIR colours into the three paramet ers returned by the fitting 
procedure: T d , B, and NdK^. We refer to lShettv etafl J2009allbh for 
a more detailed description of the issues and uncertainties in using 
simple modified blackbodies for integrated IR data. 

We fit the modified blackbody to the 100-500 jim Herschel 
bands at the 500 fim resolution. We do not include the PACS 70 /jm 
data due to the lower S/N across the M31 image, and as we expect 
the emission at these wavelengths to have a significant contribution 
from stochastically heated dust grains. We limit our fitting to all 
pixels with S IN > 5 in all 5 bands, which in practice is limited pre- 
dominantly by the PACS lOOyurn band. To fit the simple modified 
blackbody to the data, we assume a uniform, bounded prior grid 
for all three parameters, compute the x 1 goodness of fit for every 
model parameter set j, 



z 



(2) 



where 5 v , b 8 and <x,, are the observed surface brightness and ob- 
servational uncertainty for the band v, respectively, and S v j is the 
model flux for parameter set, j, determined from Equation Q] The 
bounds for dust temperature and emissivity were 5 < T A < 50 and 
0.5 < B < 3.5, respectively. We then determine the probability for 
each parameter set assuming a gaussian distribution, exp(-^/2), 
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Figure 5. Radial variation of mean temperature in M3 1 . The temperature 
is determined by fitting the median flux in all bands in radial bins of 230pc 
width using elliptical annuli. The error bars show the 16-84 percentiles of 
marginalised PDFs. 



and marginalise over the other parameters to determine the prob- 
ability distribution function (PDF) for each parameter, e.g. T d . An 
examination of the PDF for T d in several pixels shows a symmetric 
Gaussian distribution (as can be seen in the symmetric uncertain- 
ties in Figure [5), providing justification for our assumed Gaussian 
uncertainties for the parameters. 

4.2 The dust temperature across M31 

We performed this fit for every 14" pixel in the M31 IR data that 
laid above our 5 IN cut, with the resulting image for the dust tem- 
perature shown in Figure|4] Pixels which fall below our S /N cut are 
shown as white, with the dust temperature (K) shown by the colour, 
as indicated by the colour-bar at the bottom of the image. Note 
that within 2 kpc, almost all pixels have sufficient S/N. The median 
temperature across the image is 17K, close to MIPS based disk 
tempe rature of 18 ± 1 K determined by Tabatabaei & Berk huiisenl 
(2010). The uncertainties (half the 16-84 percentiles, or 1 <x if the 
uncertainties are Gaussian) range from 5 to 15%, with the largest 
uncertainties in regions at the lowest signal to noise and the central 
kiloparsec. The median uncertainty is ~6% (~ IK) dominated by 
the emission in t he dusty, star-form ing ring seen in GALEX and 
FIR images (e.g.. iThilker et al .1120051) 

Two things stand out in this temperature map: the almost con- 
stant temperature in the disk of M31, with T d = 17 ± IK, and the 
strong temperature increase in the central ~2 kpc. The temperature 
in the lOkpc ring of M31 is clearly not exactly constant, with warm 
spots occurring throughout the ring. These warm spots correspond 
with the locations of Hn regio ns in the disk as dete r mined by Ha 
image s (c.f. the Ha maps of iDevereux et al.|[l994l : lAzimlu et al.l 
l20ll . 

4.3 The radial dust temperature 

The central dust temperature rise is more clearly quantified by tak- 
ing the median SED in radial bins. To do so we determined the 
median flux in all bands in radial bins of 230 pc (using elliptical 
annuli with a P.A.= 37.7° and axis ratio of 0.26), and determined 
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Figure 4. Effective dust temperature map of M3 1 , determined from a single modified blackbody fit to each matched pixel of the PACS 1 00, 1 60//m and SPIRE 
250, 350, and 500 fim bands, with all bands convolved to the SPIRE 500 fim PSF and resolution (as described in sections l4~Tl and l4~2l . The median temperature 
is 17K, and is dominated by the emission in the ring, while the highest temperature is reached in the nuclear region, as can be seen by the clear gradient in the 
figure. The median uncertainty is ~ IK or ~6%. 



the PDFs for all three parameters for each radial bin. The result- 
ing temperature variation with radius is shown in Figure [3] with 
the 16-84 percentiles marked by the error bars. This plot clearly 
shows the marked increase in dust temperature in the inner 2kpc 
from the median disk dust temperature of T d = 16 - 17 K. The peak 
in dust temperature at the centre of ~ 35 K is similar to t he dust 
temperature of T d = 33 K determined bv lSoifer et al, for the 

central 4', using IRA S data and assum ing a power law emissivity 
of p = 2. Similarlv. lHabingetalJ d 19841) found 34K (using j3 = 1) 
for approximately the same central region and same data. As these 
works use apertures of 4' diameter, to define the centre, these tem- 
peratures represent an average of the inner ~ 0.35 - 0.45 kpc and 
thus are remarkably close to our determined value. As these fits are 
based on only the IRAS 60 and 100 jim bands, and thus biased to 
the warmer dust temperatures due to the shorter wavelengths, the 
similarity in temperatures suggest that a single, warm component 
of dust dominates the IR SED in the centre. This is supported by 
the similarity of the PACS 70 jim flux and the model flux in the 
central kpc, even though not actually used to determine the fit. 

In addition to giving access to the full IR SED, and thus a bet- 
ter measure of dust temperature, the higher resolution of Herschel 
enables us to see the steep gradient in temperature from the disk to 
the centre. This difference bet ween the bulge and disk du st temper- 
ature in M3 1 agrees with what Engel bracht et al.l d2010h found fo r 
the KINGFISH sample of nearby galaxies jKennicutt etai]|201lh . 
lEngelbracht et al.l d2010h found that the ratio of the central to disk 
dust temperature was greater than 1 across their sample, and in- 
creased with increasingly earlier types (see e.g. their Figure 3). 

The emissivity,/?, determined from these fits has a large uncer- 
tainty across M31, and is consistent with /? = 2 at all radii. Given 
this, we created a higher resolution temperature map using only 
the 100, 160, and 250 yum bands with the assumption of a constant 
emissivity of fj = 2. The resulting image of the central region is 
shown in Figure [6] where the central kiloparsec region is marked 
as a black circle and the median uncertainty is 0.75 K. Interestingly, 
while more structure in the temperature distribution in the centre is 



25.0' 



20.0' 



10.0 




30.0s 43m00.0s 30.05 0h42m00.0s 
RA (J2000) 

Figure 6. Effective dust temperature map of the central region of M3 ^de- 
termined from a single modified blackbody fit to each matched 6" pixel of 
the convolved PACS 100, 160 and SPIRE 250 fim bands, assuming a con- 
stant emissivity of /? = 2.0. The black circle outlines the central kiloparsec 
region. 



seen due to the higher resolution, the overall temperature gradient 
in the centre still dominates the image. The steeper gradient in dust 
temperature along the minor axis arises from the mixture of disk 
dust emission with that from the bulge, leading to an overall lower 
average dust temperature in these regions. The fine scale structure 
is likely due to the distribution of the diffuse gas and dust as seen 
in Figure|2] 



4.4 The heating mechanism of the bulge dust 

The clear temperature increase in the central regions of M3 1 as seen 
in Figure [5] is interesting as it closely corresponds with the stellar 
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light distribution as traced by the near-IR emission, suggesting tha t 
the heating of the dust and bulge are linked. As iLi et al.l (2009) 
demonstrated (particularly their Figure 3c), both the Ho- and dif- 
fuse X-ray emission also increase toward the centre in the same 
manner as the stellar light. However, with the resolution offered by 
Herschel we can begin to understand the heating of the dust in the 
bulge of M31, and by proxy in the bulges of early type spirals and 
in early-type galaxies. 



4.4.1 Potential heating sources 

In principle, one possibility for heating the dust is a low level of star 
formation, buried amongst the dust and bulge stars. However, high 
resolution observations by HS T show that no OB stars are appar- 
ent thro ughout the bulge (e.g. iBrown et al.lll998l : iRosenfield et al.l 
lin prep.1) . and no regions of high eno ugh attenuation to bur y such 
potential star formation are observed (Melchior et al.lcO OO). Also, 
the observed emission line ratios in the ionized gas (e.g. [N ii]/Hq') 
argue against young stars being the dom inant heating mechanism 
dCiardullo et al.ll 19881; ISaglia et all2010t) . 

As the super m assive black ho le at the centre of M3 1 (dubbed 
M31*) is quiescent dLi et alj|201ll) . heating by an AGN radiation 
field can be discounted. In addition, the radial profile of dust tem- 
perature is incompatible with a dominant central source heating, as 
shown in the following section. 

Collisional h eating of the dust by the hot 10 6 K X-ray gas at 
the centre of M31 ( Li & Wang 2007) is another possibility (see, e.g. 
1^1199 ll : lNataie" et alj2 010). However, the association of the dust 
with the Ho- emission and the weak detection of CO in regions of 
higher attenuation (discussed in section [3~2l indicate that the dust is 
in a denser medium, and not predominantly associated with the X- 
ray gas itself. Likewise, the total energy in the diffuse X-ray gas is 
insufficient to fully explain the Ho - emission, and thus is unlikely to 
heat the dust either dLi et al.l2009h . Thus, while this may contribute 
to the emission, it is insufficient to explain the majority of the hot 
dust emission. 



4.4.2 Bulge star heating 

This leaves the high stellar density of old stars at the centre as the 
dominant heating mec hanism of the dust in the bulge of M31, as 
originally suggested bvHab ing et alJdl984T) . By comparing the stel- 
lar light distribution for the bulge against the distribution of tem- 
peratures in the central kiloparsecs, we can now demonstrate the 
link between the heating of dust by the bulge and the dust emis- 
sion. In the case of a simple steady-state temperature, grain cooling 
balances the dust heating, an d the grain te mperature is then pro- 
portional to (as shown in, e.g. lDrainefeoi ll particularly Equations 
24.17 and 24.18); 



'hvo 



60T(4 + B){(4 + B)Q cr 



1/C4+/3) 



u: 



1/I4+/J) 



(3) 



where Q and (2^)* are the dust absorption cross-section at 
the reference frequenc y vp and aver aged across the heating spec- 
trum (Equation 24.2 in iDrainellioT 1), respectively, T and £ are the 
Gamma and Riemann zeta functions and U, is the heating stellar ra- 
diation field density. In the simple assumption where the emissivity 
slope is constant at B = 2, this reduces to 



\1607T 2 ^Cr 



(<2ab s >*,0V§{/,)' 



(4) 



where (2abs)*,o ls the spectrally-averaged dust absorption cross- 
section normalised at vo- Using a reference wavelength of A = 
100 fim and the unattenuated spectral energy distribution deter- 
mined in section [331 (<2 ;l bs)*,o ~ 197 for a Milky Way dust model 
with Ry = 3.10 (at A = 250 yum, (<2 a bs)*,o is approximately (2.5) 2 
larger: |Prainel2003h . 

We can estimate the stellar radiation field heating the dust with 
a few simple assumptions: 1) the dust is in a thin disk at the cen- 
tre of M3 1 and 2) is optically thin to the dominant heating radia- 
tion, and 3) the bulge stellar mass is spherically distributed and 4) 
has constant colours and stellar population. The first assumption is 
based on the findings given in section [3721 where the correlation of 
ionized gas and dust suggest a thin disk geometry. The measured 
opacit y in the bulge of M3 1 is low, as shown by iMelchior et al.l 
(2000), as is the average dust column of ~ 4.5 x 10 4 M o kpc 2 
(see 93. 3K thus the assumption of optically thin dust is not un- 
reasonable. As for the third assumption, the bulge of M31 has al- 
ready been shown to be non- symmetric, displaying a boxy struc- 
ture, and possibly even differen t RA. and inclination angle to the 
disk (e.g. ICourteau et alj|201lh . However, a simple spherical ap- 
proximation for the bulge with a constant mass-to-light ratio can ac- 
tu ally reasonably well reproduce the stellar light profile, as shown 



tually reasonably v 
bv lGeehanetalj iB 



2006). The NIR colours are close to constant as 



shown in Figure 1 in Courteau et al.l d201 ll) . but th e UV colour, 
howe ver, does show a gradient (see, e.g. Figure 3 in lThilker et al.l 
2005). While some of this may arise from a gradient in attenuation, 
there e xists also radial gradients in the UV-emitting stellar popu- 
lations dRosenfield et alJlln prepj) that may affect our assumptions 
(though the follow section argues against this). 

As shown in Appendix [A] given these assumpti ons, using the 
model for the bulge profile from Gee han et al.l d2006l) and the inte- 
grated unattenuated stellar luminosity from section 13.31 the bulge 
interstellar radiation field heating the dust is 



U,(r) ■ 



3.66 x 10 



Jo 



dx erg cm 3 , 



(5) 



(r/kpc) Jo (1+x) 3 

at a radius r (in kiloparsecs), and a bulge scale radius of r b = 
0.61 kpc. 

From Equations [4] and [5] we can then determine the expected 
dust temperature (given the above assumptions) as a function of 
spherical radius. For an infinitely thin disk, located in the plane of 
M31, we can equate the dust temperature in spherical coordinates 
with that in cylindrical (i.e. r = R, z = 0), giving 



T<i,v,(R) = 1750 [U,(R)]~ 6 K, 



(6) 



where the factor 1750 arises from using the Milky Way dust model 
described above for (l2abs)*o, and U, is in ergcirT 3 as in Equation 
® 

In Figure [7] we show the distribution of the 6" pixels from 
the dust map in Figure [6] in terms of the dust temperatures with 
distance from the centre, based on our simple circular radius. Over- 
laid on this distribution are two curves. The expected dust heat- 
ing distribution, T iVt (R), as given by Equation [6] is shown by the 
dotted curve. As is clear from the figure, this is in excess of the 
observed temperature in M31's bulge. The solid curve shows the 
same r dUt (K), however multiplied by a factor 0.8. This brings it 
in agreement with the observed temperature distribution except at 
radii larger than ~ 1.2 kpc. The first thing to take from these two 
curves is that the bulge provides a more than sufficient interstellar 

4 Available fromhttp: //www. astro .princeton. edu/~draine/dust/dustmix. html 
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Figure 7. The distribution of dust temperatures within the inner 2 kpc 
(530") of M31. The colours show the number density of the 6" pixels from 
Figure[6]in terms of (T$) and circular radius, as labelled by the colour bar to 
the right. Overlaid on this are two curves showing expected heating as based 
on the determined bulge ISRF. The dotted and solid curves show, respec- 
tively, the expected dust temperature distribution, 7^ u, (R), from Equation 
|6]and the same 7j u„(S) scaled by 0.8. 



radiation field (ISRF b ) to heat the dust, with expected temperatures 
in excess of that observed. However, importantly, the observed gra- 
dient in temperatures matches that expected from heating the bulge. 
This leaves the question, why is the expected temperature distribu- 
tion in excess of that observed? 

There are several possibilities given our simple assumptions. 
One possible explanation is suggested by the flattening of the r dlist 
pixel distribution in Figure [7j at radii > 1.3 kpc. This flattening is 
most likely due to the increasing disk contribution to the IR at these 
radii, which acts to lower the mean line-of-sight dust temperature. 
This is approximately the radius where the stellar mass distribu- 
tion , and therefore light distribution, becomes disk dominated (see, 
e.g. iGeehan et al1l2006l : ICourteau et alll20l ]]), and thus disk dust 
is likely also dominating the determined dust temperature. As dis- 
cussed in section 1431 and shown in Figure [4] the mean dust tem- 
perature across the disk is ~ 17 K, which is also where the pixel 
distribution flattens to in Figure [7] Figure [6] also gives some indi- 
cation of the disk contribution to the bulge temperature, as there 
is clearly a flatter gradient of temperatures along the major axis as 
compared to the minor axis, which should have a greater contribu- 
tion of disk dust due to the high inclination of M3 1 . However, an 
examination of the observed IR SED within the central kiloparsec 
indicates that it is well described by a single black-body, as men- 
tioned in section |4~3l As a confirmation we also fitted the pixel- 
SEDs in the central region with a combination of two modified- 
blackbodies. For both blackbodies, we assumed an emissivity of 
[3 = 2, and let the temperature of the first component free, and that 
of the second, cooler component constrained to 17 K. A simple fit 
found that the central kpc region is dominated by the warmer com- 
ponent, whose temperature distribution is almost indistinguishable 
from FigureUJ Thus, while it must occur at radii > 1 kpc, the contri- 
bution of the disk cannot explain the offset between the theoretical 
curve and the observational-based temperature. 

Another possibility is that the dust in the bulge is not Milky 
Way like, and our assumed (<2abs}*,o is incorrect. For example, if 
we assume a Milky Way model of dust, but with R\ = 5.5 (using 



the Jweingarfner & Drainel2 001) model opacity), <2abs>*.o reduces 
to ~ 155 (normalised at Ao = 100 //m, as above). Such a possibility 
is reasonable, given that the extreme environs of the bulge may act 
to alter the dust size distribution, destroying small grains leading 
to a flatter opacity. However, given the index of 1 /6 in Eqn. [4] the 
offset in temperature requires a significant change in dust opacity 
if this alone causes this offset. Similarly, our estima te of U, is also 
likely incorrect given that the bulge is not spherical dCourteau et al.l 
1201 lh . and the dust is not in a perfectly thin disk. However, the 
index of 1 /6 again requires our U, to be over estimated by approx- 
imately an order of magnitude if it alone is incorrect. 

More likely, it is a combination of these, plus the possibility of 
some self-shielding by dust (which acts to reduce the U, seen by the 
dust), which lead to the difference between the theoretical tempera- 
ture distribution and that determined from the modified blackbody 
fit to the data. Yet it should still be noted, that these are all needed 
to reduce the theoretical heating to that observed. The bulge stel- 
lar radiation field provides more than sufficient energy to heat the 
warm dust at the centre of M31. 

As a final check that the bulge stars are dominating the heat- 
ing, we can apply the same methodology and reasoning as used 
for Eqn. [4] but assuming that the radiation field heating the dust 
arises from a central point source (i.e. AGN or nuclear star cluster), 
meaning that U oc r 2 . However, when T iVt (R) is determined for 
this radiation field (similar to Eqn. |6}, the slope does not show the 
same form as in Figure UJ demonstrating that the heating radiation 
field must be extended. Even allowing for a greater contribution of 
the disk emission to approximate the observed slope, to match the 
point source T i Ut (R) to the measured T A requires a central luminos- 
ity that would be obvious in the optical (or X-ray) emission, whic h 
is discounted by observations as already discussed dLi et al.|[2009h - 



4.5 The heating of dust by old stars 

While Figure |7J clearly links the bulge with the dust emission, it 
does not reveal what exactly is heating the dust. Typically, due to 
the steep wavelength dependence of the dust opacity, the dust IR 
emission has been associated with UV light and hence star for- 
mation. Typically, diffuse dust is considered to be heated by a lo- 
cal interstellar rad iation field (ISRF) or scaling thereof (see e.g. 
Draine &Lill2007t). As this I SRF is typically based on the local 
ISRF of lMathisetalJ ll983l), the UV light do minates the heat- 
ing of dust (e.g. Figure 4 in Mathi s et allll983l) . and as the local 
UV light is dominated by B stars, the diffuse dust emission in- 
directly traces B stars and thu s star formation on longer (~ 100 
Myr) timescales. As shown bv lThilkeretal] J2005I) . the UV light 
strongly peaks at the centre of M31. However this UV light in 
the bulge has been shown to be not associated with st ar-formation, 
but rather arises from extreme h orizontal branch stars (brown et al.l 
ll998l;lRosenfield et alllin prep.l) . 

However, while these hot, low-mass stars may be the source 
of the UV upturn and possibly the Ho- emission, they cannot be 
the dominant heati ng source of t he dus t. Based on the IRAS and 
ANS observations, lHabing et al] d 19841) realised that the UV light 
would not be sufficient to heat the dust, and the heating would be 
dominated by the A > 300 nm light from evolved stars. In Figure|8] 
we show three panels that illustrate which radiation heats the dust. 
In the top panel we take the "unattenuated SED" from Figure [3] 
which is likely an overestimate at A < 0.25pm (see Section[3}. The 
middle panel shows the total dust absorption cross-s ection per gram 
of dust for the MW dust model with R w = 3. 1 from Draine ( 2003), 
as also used in the previous section. By multiplying these together 
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Figure 8. The energy absorbed by diffuse dust in the bulge of M3 1 . The top 
panel shows the "unat tenuated" bulge UV-NIR SED from the MAGPHYS 
Ida Cunha et alj2008t) model fit to the integrated SED of the central kilopar- 
secs (blue curve in Figure |3}. The middle panel shows the total absorption 
cross-section per gram of dust for the MW dust model with i?v =3.1 from 
iDraind j20O3h . The bottom panel shows the product of these, revealing the 
wavelength distribution of the energy absorbed by dust. The conclusion is 
that the optical stellar light is the main contributor to the dust heating in the 
bulge. 



we obtain the amount of energy absorbed per H atom as a function 
of wavelength (bottom panel). This Figure reveals that most of the 
energy that heats the dust arises from optical photons, not the UV. 
As the Q a bs( v ) slope is « 1, the UV-optical slope needs to be < -1 
(in frequency) for optical photons to dominate the heating of dust. 

As mentioned before, the bottom panel of Figure[8]most likely 
overemphasises the role of UV photons in the heating due to the 
small UV excess in the model stellar SED (blue curve in Figure[3]l. 
The intrinsic UV emission produced by the evolved stars in the cen- 
tre of M3 1 (the so-called UV- upturn observed in early-ty pe galax- 
ies) is well-reproduced by the iBruzual & Charlotl d2003l) code, as 
shown by a comparison of models with the obs erved UV colours o f 
early-type galaxies (see specifically Figure 5 in lDonas et all2007r) . 
Even if the model UV is incorrect, the intrinsic UV flux would still 
need to increase by ~ 0.4 dex (or 0.7 dex from the observed value) 
to be commensurate with the energy deposited in the optical. This 
lack of UV heating is also cl ear when the top pa nel of Figure [8] 
is compared with the ISRF of iMathis et alj J 1983!) (e.g. their Fig- 
ure 1). When normalised to the peak at ~ 1 yum (in both radiation 
fields), the UV in the IRSF is approximately an order of magnitude 
greater than that determined in M3 1 . 

Furthermore, the hot environment of the circumnuclear region 
of M31 may destroy the smallest dust grains as discussed in the 
previous section. The destruction of small grains flattens the slope 
of /Cabs(v), leading to yet a greater dominance of optical photons 
in the dust heating. Thus the bulge of M3 1 is a clear case where 
the heating of dust is dominated by the light from old stars. This 
is in stark contrast to typical star-forming galaxies, whe re the UV 
from young stars dominates the heating of dust (e.g. lLaw et all 



1201 ll) . iMontalto et all d2009h had already found with GALEX- 
Spitzer data that the dust in the Andromeda galaxy a s a whole ap- 
ears to be heated mainly by stars a few Gyrs old. lHabing et al.l 
1984) already inferred that the bulge of M31 takes this further 
with heating likely dominated by evolved stars. However, the mean 
stellar age of the bulge of M3 1 is > 6Gyr, with some e s timates 
placing the dominant ag e > 10 Gyr (Davidge e t alj|2005l : iBrownl 
2009; ISaglia et"al]|2010l) . demonstrating that even the oldest stars 
are able to heat dust. This is a clear cautionary note against using 
warm IR emission as a direct tracer of star formation, or obscured 
UV emission. 



5 DISCUSSION 

As the nearest massive galaxy, the Andromeda galaxy allows us 
to connect the small scale physics with the integrated properties 
of galaxies. The central kiloparsec of M3 1 actually matches very 
well many of the observed properties of nearby early-type galax- 
ies. Within a 1 kiloparsec (4.4') circular aperture the estimated 
stellar mass is ~ 10 10 M o , with a ver y old, red stellar l ight, esti- 
mated to be greater than 10 Gyrs old (Saglia et al. 2010). The op- 
tical colours of the bulge, e.g. NUV-r x 5.0, pla ce it well in the 
realm of gas poor early-type galaxies (see, e .g. lOke & Sa ndage 
1 19681 ; ISaintonge et alj201 ll i lsmith et alj201 lh . 

There is a low level of dust attenuation across the centre, visi- 
ble in both the difference between the observed and modelled in- 
trinsic SEDs in Figure [3] and in the Ab map (their Figure 1) of 
llVIelchior et alj < l2000h . This low level of attenuation is a result of 
the low dust column across the centre and connected to the rel- 
atively weak IR emission. The tot al dust mass within this aper- 
ture from the lda Cunha et al.ld2008l) MAGPHYS model is 10 5 2 M G , 
contributing only 0.5% of the total dust mass to M31. However, 
the bulge contributes ~ 5% of the total IR luminosity due to 
the relatively warm, blue FIR emission. This dust mass results 
in a very low M dl , st /M t for the bulge, with ~ 10~ 47 . This actu- 
ally places the bulge of M31 in a similar regime as the sample 
of earl y-type galaxies exp lo red in the Hersche l reference survey 
(HRS, [Boselli et alj20ld) in l Smith et ail feoill) . falling somewhat 
in between the SO and E galaxies (see their Figure 8). Similarly, 
the temperature increase from disk to bulge in M31 follows the 
tre nd of having warmer dus t in earlier Hubble-types, as observed 
bv lEngelbracht et all d2010h in the KINGFISH sample. This trend 
was also seen by Smith et al. I J20T ]]) in the HRS, who found warmer 
dust in E galaxies than SO, and an overall warm dust temperature 
for the sample (mean Ta = 24K). 

The M gils /M t ratio is also very low in the bulge of M31, with 
little Hi and CO detected in the centre. However, ionized gas is 
seen in Ha and the FIR emission is spatially well-correlated with 
this gas showing a similar lower inclination, barred-spiral pattern 
as visible in Figure [2] with the MIR e mission also following the 
diffuse Ha morphology dLi et alj|2009l) . Given this correlation, the 
Ha- emission allows the structure of the dust distribution to be de- 
termined (e.g. Ilacobvetalj 19851) . and indicates that the ISM in the 
centre, including the dust, is in a thin, spiral disk. The origin of 
this gas and dust in the cen tre of M31 is s till not known, and be- 
yond the scope of this work. lLTet al.l 12009) suggest that the stellar 
ejecta are more than sufficient to replenish the observed hot gas in 
M31's centre, and suggest the rest of this matter comes out as a hot, 
X-ray flow. However, the spiral pattern observed in the bulge does 
appear to link up with the emission in the disk (Figure [T), even if it 
appears to be at a different inclination to the disk. 
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Given the above similarities using the bulge of M3 1 as a re- 
solved representative of the spheroids of other early-type galaxies 
is a reasonable assumption. Based on this and our resolved study of 
dust heating in the bulge we can extrapolate to infer that a similar 
heating mechanism for the dust occurs in th e majority of earl y-type 
galaxies in which dust is detected (i.e., the lSmifh et alj ( 1201 lh sam- 
ple). While many early-type galaxies do have some AGN activity 
in the centre, and several, generall y with higher gas ma sses, have 
observed active star formation (e.g. ICrocker et al.lbOllb . the com- 
bination of diffuse, optically thin dust and a strong, diffuse radiation 
field due to the high stellar density is able to heat dust at the cen- 
tres of these early-type galaxies to a significantly warmer level than 
that observed in most stellar disks. This increased dust temperature 
makes the generally dust poor early-types still visible in the recent 
and ongoing FIR surveys with Herschel. 

However, the possible contribution of disk dust to the observed 
IR emission, as derived from the dust temperature offset in Figure|7] 
suggests a final cautionary note in this extrapolation. The bulge of 
M3 1 lies in a far different environment to typical early-type galax- 
ies, which affects both the observed emission, and the evolution of 
the ISM in this early-type spheroid. 



6 SUMMARY 

When observed in the far-infrared, the bulge of the Andromeda 
galaxy (M3 1 , NGC 224) stands out as a region of luminous blue 
emission (i.e. 70 bright), surrounded by ring-like red emission 
from the disk. This corresponds to a peak in the mean dust temper- 
ature, Td (Figure[4}. Across the disk of M31, the mean dust temper- 
ature is reasonably constant with Tj = 17 ± IK. However, within 
the central kiloparsec, the temperature of the dust rapidly increases, 
reaching Td ~ 35K at the centre, as seen in Figure[5] 

The start of this upturn is also where the bulge begins to dom- 
inate the stellar radiation field, clearly indicating that the heating 
mechanism of this warm dust and the bulge are associated. Heating 
by either an AGN or gas-grain collisions is discounted by X-ray 
observations, leaving stellar heating as the likely source. The radial 
profile of this steep increase in temperature in the inner kiloparsecs 
corresponds well with the temperature slope expected from heat- 
ing by the diffuse radiation field arising from the bulge stars, sug- 
gesting direct link. The theoretical dust temperature based on the 
diffuse radiation field and our assumptions is actually in excess of 
that observed, providing more than sufficient heating. This excess 
suggests the not unreasonable possibilities of different dust proper- 
ties in the bulge from the standard diffuse dust in the Milky Way 
and that the dust is likely to be in a clumpy distribution with some 
self-shielding occurring. 

However, even though the bulge of M31 is observed in the 
UV, no young stars are seen in HST ob servations or spectra 
dRosenfield et aljin prep.l : ISaglia et al.l2010 ). leaving old (> 6 Gyr) 
stars as the dominant heating mechanism. Furthermore, by taking 
the observed SED of the bulge and standard dust opacity it is clear 
that it is the red optical light of these old stars that dominates the 
heating of the warm dust in the centre, as originally suggested by 
Habin getail jl984h with the IRAS observations. 

Together this demonstrates the possibility of heating dust by 
stars > 10 billion years old. This is one of the clearest demonstra- 
tions of how IR emission does not always correlate with star forma- 
tion, a common assumption, but rather is dependent upon the dis- 
tribution of dust (and associated gas) and the total radiation field, 
which includes both young and old stars. The bulge of M3 1 also 



demonstrates one of the plausible mechanisms for heating the rela- 
tively warm dust observed in the centres of early-type spirals, and 
early-type galaxies. 
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d2006l) and use a iHernquistl ( fl99(jl) mass profile due to its simple 
spherical geometry and assume a constant mass-to-light ratio, giv- 
ing; 

L b 1 



v*(r) 



(Al) 



iGeehan et all j2006h found the bulge radius to be 0.61 kpc (their 
Equation 1 and Table 2), and we used the determined unattenuated 
stellar luminosity within 1 kpc of L t =~ 10" L Q (from section [3~3l 
to normalise L b to 10' 032 L o - 

The interstellar radiation field in the bulge (ISRF b ) is then sim- 
ply the luminosity density (in e.g. L Q pc 3 ) convolved with spherical 
dilution, 1/r 2 ; 

1 



U,{r) 



Anr 2 c 



v*(f"o)- 



1 



-dV, 



(A2) 



An(r -f) 2 c 

where the factor 1 jc is introduced to convert the ISRFb to an energy 
density (i.e. erg cm -3 ). As both the luminosity density and dilution 
are purely radial functions, the volume integral in Equation lA2l re- 
duces to (in spherical coordinates); 



U,(r) 



r y»(ro) r 2 " r 

Jo 4;rc Jo Jo 



■ 2rrn cos d> + r 2 



dipdddro 



2rc V |r — rb| 

Substitution of Eqn.[AT]into this and condensing then gives 



(A3) 



U«(r) = dx, 

Anr b c r Jo (1 + x) 3 



(A4) 



where x has been substituted for r /r b . The factor in front of 1/r 
and the integral equals 3.66 x 10~" ergcirT 3 for our given L b and 
r b , with r in kiloparsecs. While the integral in Eqn. lA4l can be ana- 
lytically determined, we use a numerical integration for simplicity 
within this work. 



APPENDIX A: THE BULGE RADIATION FIELD 

For the bulge radiation field we first need the stellar l uminosity dis- 
tribution (which we label v+(r)). For this we follow IGeehan et all 
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